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We generalize the topological response theory of three-dimensional topological insulators (TI) 
to metallic systems - specifically, doped TI with finite bulk carrier density and a time-reversal 
symmetry breaking field near the surface. We show that there is an inhomogeneity-induced Berry 
phase contribution to the surface Hall conductivity that is completely determined by the occupied 
states and is independent of other details such as band dispersion and impurities. In the limit of 
zero bulk carrier density, this intrinsic surface Hall conductivity reduces to the half-integer quantized 
surface Hall conductivity of TI. Based on our theory we predict the behavior of the surface Hall 
conductivity for a doped topological insulator with a top gate, which can be directly compared with 
experiments. 



Topological insulators (TI) are states of matter with 
a bulk energy gap and topologically protected gapless 
edge states. Their most striking properties are their 
topological responses, such as the quantized Hall conduc- 
tance of integer quantum Hall (QH) states [1]. Recently, 
a wide class of time-reversal (TR) symmetry protected 
topological insulators has been discovered [2-5]. The 
three-dimensional ones are crystalline materials with a 
bulk energy gap and massless linearly dispersing Dirac- 
like surface states [6-8] . The topological responses of 3D 
TI can be observed when the surface states are gapped 
by magnetism [1] or superconductivity [10]. When the 
surface states are gapped by magnetism, the TI obtains 
a topological electromagnetic response described by the 
effective action = J d'^xdt^ ■ B with a = /he 
the fine structure constant, and = tt modulo 27r [1]. 
This is a topological magneto-electric effect and can be 
observed as a surface QH effect with a half integral quan- 
tized Hall conductivity axy = ~ (^^^)7r' ^^^^ ^ 
an integer. For non-interacting topological insulators 
can be calculated explicitly as the Chern- Simons invari- 
ant of the geometrical gauge field defined in momentum 
space [1, 11]. The topological response theory of 3D TI 
can be obtained by a dimensional reduction from the 4D 
generalization of the quantum Hall effect [12]. 

With the recent systematic understanding of topolog- 
ical insulators, a natural question is whether only insu- 
lators can be topological: are there topological phenom- 
ena in gapless systems such as a Fermi liquid? In this 
paper, we develop a theory to characterize the intrinsic 
response properties in a generic 3D Fermi liquid. Al- 
though in a metallic system the response properties such 
as Hall conductivity are not quantized, there can still be 
intrinsic "topological" contributions that are determined 
by the Berry phase gauge field in momentum space and 
are independent of non-topological details such as energy 
dispersion away from the Fermi surfaces and impurities. 

An interesting example of a 3D Fermi liquid with non- 
trivial topological response is a doped Tl-a system with 
the band structure of a TI but with a finite bulk car- 




FIG. 1: Schematic picture of a doped TI with a ferromagnetic 
layer on the surface. The spheres indicate the Fermi surface 
size which is position dependent. Near the surface, the Fermi 
level lies in the gap. 



rier density. Here, we study the surface Hall conduc- 
tivity of a doped TI with TR breaking field applied at 
the surface, as illustrated in Fig. 1. (Some other be- 
havior of doped TI inherited from the nearby TI has 
been studied recently [13, 14]). Unlike a TI, the sur- 
face Hall conductivity of a doped TI is not quantized, 
and is in general dependent on surface conditions. How- 
ever, in the smooth boundary limit, where the changes of 
the Fermi level and the band structure near the bound- 
ary are smooth on the scale of the mean free path, we 
show that the surface Hall conductivity is determined 
by a topological invariant — the Chern- Simons invariant 
of the geometrical gauge field integrated over the two- 
dimensional Fermi surface and the real space direction 
perpendicular to the boundary. We obtain our result by 
generalizing the dimensional reduction approach in Ref. 
[1] to metallic systems. Our result can be understood 
as a higher dimensional, inhomogeneity-induced gener- 
alization of the intrinsic anomalous Hall effect [15-18] 
and its interpretation as a topological property of Fermi 
liquids [19]. Our results for doped TI are relevant to the 
current experiments on TI materials since all the known 
3D TI materials, such as Bii-a^Sb^^ alloy [20, 21], Bi2Se3 
and Bi2 Tea [22-24], have finite residual carrier density, 
which has become a major experimental challenge for 
the observation of topological effects in TI. We also nu- 
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merically studied the surface Hall conductivity of doped 
Bi2Se3 with a simple surface T-breaking field using the 
four band effective model [22, 25]; our results can be ver- 
ified by zero-field transverse conductance measurements 
of Bi2Se3 with a magnetic layer deposited on the surface 
and voltage applied at a top gate. 

Consider a Hamiltonian Hsd[A^^ ^(^)] foi" ^ 3+lD non- 
interacting fermion system that is coupled to the electro- 
magnetic gauge field A and that depends on a spatially 
varying parameter 0{r). For example, in a doped TI with 
a boundary, 0{r) may describe both the spatial variation 
of the Fermi energy and the change of band structure 
across the boundary. Without losing generality, we can 
define ^ = and tt to describe the vacuum and the bulk 
of the doped TI, respectively. To study the surface of a 
doped TI, we take 0{r) = 0{z) to depend on one direc- 
tion of space, with 6{z ^ — oc) = tt deep in the bulk, 
while ^(0) = at the surface. When varies smoothly, 
at each point in space the system without external mag- 
netic field is locally described by a Bloch Hamiltonian 
H = ^^c\.h[k^6{z)\c\f^^ where k = (kx^ky^kz). For N 
bands, Ck is a A/'-dimensional vector and h[k,0{z)] is 
di N X N Hermitian matrix. In the physical system 
is defined in the region [0,7r]. It is convenient to ex- 
tend the definition of h[k^O] to G [0, 27r] by setting 
h[k,0] = rt/i[-ifc,27r - 0]T for tt < 6> < 27r. Here T is 
the time-reversal operator, and h[k,0] is continuous in 
since it is TR invariant at ^ = 0, tt. 

Now we study the electromagnetic response of the sys- 
tem, which is determined by the effective action 

^iS,ii[A^,e] ^ J f dt[Y:iCiidtCi-H3D[A,e]] ^ (^-^^j 

The DC Hall conductivity is determined by the quadratic 
term in the action: 

SfiM^,0] = ^J d^xdte^'"^-a^[e]A,d„Ar. (2) 

with the Greek letters /i, z/... = 0,1,2,3 labeling the 
time and space coordinates. The spatial components 
cr^, i = 1,2,3 of the four vector are the Hall con- 
ductivity in yz, zx and xy plane, respectively; that 
is, (Ji = e^^^Gjk^ with cjjk the conductivity tensor. In 
general, (Ji = <Ji[6>(r)] is a functional of the inhomoge- 
neous ^(r). For a smooth function ^(r), one can expand 
e{r) = 6>o + r • VO{r)\r=Q with Oq = 0{r = 0). In this 
expansion cr^ can be expanded to [31] 

ai[0{r)] ^ Cu{Oo) + ^Gs{Oo)diO. (3) 

The coefficients in the two terms are given by two- 
point and three-point correlation functions, respectively. 
Using the time-ordered Green's function G{k^uj) = 



[ct; — e , the coefficients can be expressed as: 

Cu{0,) = \j ^^e^3kir^T[GdjG-'GduG-'GdiG-'] 

(4) 

G3(^o) = -'-^e^^'' j I^Tr [Ga.G-^Ga,G-^ 

• GdkG-^GdiG-^GdeG-^] . (5) 

with dj = d/dkj and kj is the frequency and momentum 
vector with ko = u the frequency and {ki,k2,ks) = k the 
momentum. 

Before evaluating the correlation functions, let us dis- 
cuss the criteria for the validity of the above effective 
theory approach. In the case of a band insulator with a 
gap Eg^ there is a length scale ^ = hv/Eg where is a 
typical velocity scale contained in the coupling between 
conduction and valence bands. The approach above is 
valid as long as the characteristic scale of variations of 
is much longer than ^, i.e., |V^| <C 1/^. For a perfect 
metal, the gap vanishes and the Green's function in (4) 
and (5) is singular. With impurities, the Green's func- 
tion obtains a finite self-energy l/r, which corresponds 
to a length scale I = vpr, which is the mean- free path. 
The mean-free path plays the role of ^ in the insulator 
case, and the expansion above applies to the system with 
|V^| <C l/L Physically, the expansion applies because 
the conductivity at position r will be insensitive to small 
changes in that occur at r' when |r' — r| >> /. In other 
words, the intrinsic contribution to the Hall conductiv- 
ity calculated above applies in the "dirty limit" where 
the Hamiltonian is smoothly varying on the scale of the 
mean-free path, although the value of the intrinsic Hall 
conductivity is independent of the impurity strength in 
this limit. In general, the anomalous Hall conductivity of 
a metal also receives extrinsic contributions arising from 
impurity scattering, which are termed the side-jump and 
skew-scattering contributions [18]. Here we focus on the 
intrinsic contributions, which depend only on the nature 
of the Bloch states and their associated Berry phases. In 
the case of ferromagnetic metals, the intrinsic contribu- 
tion can dominate the full anomalous Hall conductivity 
in certain intermediate regimes of disorder [18]; thus we 
expect that the intrinsic contribution that we calculate 
here for the surface of a doped TI may also be the dom- 
inant contribution in certain regimes of disorder. 

Under the condition discussed above, the gradient ex- 
pansion can be done locally around each point in space, 
leading to a local Hall conductivity. For simplicity, we 
consider ^(r) = 0{z) only varying along the z direc- 
tion. The Hall conductivity in the xy plane is given by 
az{z) = Cl[0{z)] + ^[0{z)]dzO{z). As shown in Fig. 
1, we consider a single interface at 2; = between the 
vacuum at 2; > and a doped topological insulator at 
z < 0. One can always take = — tt for the bulk doped 
TI and ^ = for the vacuum, both of which are TR 
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invariant. In such a situation, the total current in the 
xy-pldLiie win be given by = / dzj^, from which we 
may obtain the two-dimensional Hall conductivity = 
J dz(Jz{z) — + <j^^'^^\ The first term = 

JdzCuiOoiz)] = ^Y.a!d^kdzF^y{z)n^{k,z) is the 
intrinsic anomalous Hall conductivity of a homogeneous 
system[15, 16], where a labels the occupied bands |ak), 
and and J^^y are the occupation number and Berry 
curvature of band a. More explicitly J^^y = dxAy —dyA^ 
with = —i{a]^\dki Modulo a quantized con- 

tribution n/27r with integer n, a^^'^^^ is determined by 
the Berry phase gauge field at the Fermi surface[19]. 
The second term al^'^'^^ = ^ J dzGs[6o{z)]dzO = 

G3{0)de is the new contribution induced by spa- 
tial inhomogeneity. With TR symmetry only broken 
at the interface, one can extend the definition of to 
[0, 27r] as mentioned above. By TR symmetry it can 
be proved that GsiO) = Gsi-O), so that dOGsiO) = 

IJ^-dOGsiO). 

Therefore, as long as the start and end points are de- 
scribed by time-reversal invariant Hamiltonians, we find 
that the integrated current response receives a contribu- 
tion from: 

^2I,;(2)^^2Z.;(2)^£!l^^ (6) 

where / = J^^ Gs{0)dO, and we have recovered the units 
e'^/h. From the definition of G3 in Eq. (5) one can see 
that plays exactly the same role as momentum and fre- 
quency. The integrand in Eq. (5) has a symmetric form 
if we view ^ as a momentum component in an "extra di- 
mension." In the case of an insulator, / is a topological 
invariant which is equal to the second Chern number as- 
sociated with the Bloch states of the 4+ ID Hamiltonian 
H[k^O]. For a TI, / will be odd, so one obtains a half- 
integral quantized transverse conductivity at the surface 
of a TI, when the surface is fully gapped [1]. 

For a doped TI, the Green's functions have poles on 
the Fermi surface, so that / is in general not quantized. 
However, / is still invariant upon arbitrary smooth defor- 
mations of the Green's function Gik^O^uo) as long as the 
deformation vanishes at the Fermi surface. Making use 
of such topological invar iance, we are able to integrate 
over frequency uj and obtain the following expression for 
/ (see Supplementary Materials for the proof): 

f = ^ [ r dOTT{TabJ'cd)e^''''. (7) 

Jbz Jo 

Here Tab is the Berry curvature associated with the bands 
that are filled, with a^b = 1,2,3,4 labeling kx^ky^kz^O. 
More explicitly J^^f = S^^f - d^Af + A]"^, 
iA^^ = {ak6\da\Pk6) with a, /3 running over filled bands. 
Formally, (7) looks the same as the second Chern number 
of the Berry phase gauge field that appears in the TI case. 




FIG. 2: Left: blue indicates schematically the band bending, 
displaying metallic bulk but gapped surface. Red indicates 
variation of 6 from to tt. Right: numerical results for a^y 
as a function of the top gate carrier density Ug = CVg/e. The 
penetration depth of the T-breaking perturbation is chosen 
to be ^ = 1 nm, the depletion length is Zd = rig /rib, the 
bulk carrier density rib = 10^^ cm~^, and the Fermi level is 
chosen to be 18 meV from the conduction band bottom [26]. 
Starting from the bottom, the different curves correspond to 
mo = 50, 100, 150, 200, 250 meV, respectively The gate 
carrier densities rig ^ 10^^ cm~^ are within experimentally 
accessible values. 

However, for metallic systems there are different numbers 
of filled bands inside and outside each Fermi surface, so 
that the definition of Tab depends on the position of the 
Fermi surface and therefore / is not a second Chern num- 
ber. For example, if there is a single Fermi surface, then 

where FFS denotes the filled Fermi sea, J^^^^ is the Berry 
curvature associated with filled bands inside the Fermi 
surface, and J''^^^ is the curvature associated with filled 
bands outside of the Fermi surface. 

Since the second Chern form is a total derivative, / 
can be expressed, modulo an integer, in terms of integrals 
over the Fermi surfaces. We find that the fractional part 
of / is: 

/modl = ^^ CcsiA""^). (8) 

where Ccs{A) = e'^^Tr[AidjAk + fAiAjAk] is the 
Chern- Simons form. The sum is over the different Fermi 
surfaces, where FSi labels the ith three-dimensional 
phase-space Fermi surface, and A^' is the Berry con- 
nection associated with the partially filled bands that 
are crossing the ith Fermi surface, i.e. that have zero 
energy for (k^O) G FSi. Eq. 8 is a consequence of a 
highly non-trivial cancellation of cross-terms associated 
with the non-Abelian CS forms, which is why only the 
partially filled bands participate. The proof of Eq. (7) 
and (8) is given in the Supplementary Materials. 

Eq. (8) is an integral over the Berry connection in- 
volving only states at the Fermi surface. This is in direct 
analogy to Haldane's result [19] that the fractional part 
of the intrinsic contribution to the anomalous Hall effect 
is due to the adiabatic phase acquired by quasiparticles 
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propagating on the Fermi surface. We see that inhomo- 
geneities can also induce Berry phase contributions to 
the anomalous Hall effect, the fractional part of which 
depends only on states at the Fermi surface. 

Note that in order for the above results to hold, we 
have assumed that varies slowly on the scale set by 
the mean free path. We must also assume that the Fermi 
surface vanishes near the surface of the material (see Figs. 
1 and 2); otherwise, the surface will cause a sharp jump 
in the Fermi level and cannot be treated as a smooth 
change. If the Fermi level already lies in a gap at the 
surface, then the sharp jump in the chemical potential 
will not have consequences for the low energy physics. 

Let us now apply our result to calculate the surface 
Hall conductivity of a realistic doped TI — doped Bi2Se3. 
From the discrepancy between bulk transport measure- 
ments like Shubnikov-de Haas oscillations and surface 
measurements like ARPES, it has been noticed that the 
energy bands of Bi2Se3 (and Bi2Te3) can bend upwards 
near the surface, so that the Fermi level lies in the en- 
ergy gap near the surface but in the conduction band 
in the bulk of the material (see Fig. 2) [26]. Such a 
situation can be further tuned using a top gate voltage. 
This can be described by our current approach since the 
gap near the surface effectively avoids the abrupt change 
in the band structure at the surface. The bulk of the 
system is TR invariant and can be described by the four- 
band effective model [2, 22] Ho[k] = e{k)I^^^ d^a\k)T'' 
with r^, a = 0, 1, 2, 3 the Dirac matrices and all the pa- 
rameters Sa\k) and e{k) given in Ref. [2] (see also the 
Supplementary Materials). 

We consider a T-breaking term induced by, e.g., mag- 
netic dopants deposited at the surface. Recently, angle- 
resolved photon emission (ARPES) experiments have ob- 
served a gap in the surface states of Mn doped Bi2Te3 
[28], with the surface Fermi level lying in the gap. This 
system may provide an experimental realization of the 
surface TR symmetry breaking, although further verifi- 
cation of the nature of the surface gap is still needed. 
This can be described by adding the following surface 
term to the Hamiltonian: 

Hi = eV{z) -h mo sm[0{z)]T^ - Mo(l + cos[(9(^)])^^ 

(9) 

where eV{z) = ~ 4^^^ ^d(l ~ ~)^ is the potential induced 
by the top gate voltage, C is the capacitance per unit area 
of the dielectric separating the top gate and the surface 
of the TI, rib is the bulk carrier density of the TI, Vg is 
the voltage applied at the top gate, Zd = CVg / erih is the 
depeletion width induced by the top gate, and e is the 
relative dielectric constant of the TI [29] (see Supplemen- 
tary Materials). Note that charge neutrality requires Zd 
to be determined by the gate voltage, which reduced the 
number of free parameters in our calculation, allowing us 
to determine the Hall conductivity from realistic material 



parameters. 

We assume that the band bending near the surface is 
entirely induced by the top gate. 0{z) = tt{1 — e~\^\^^) 
describes the surface TR symmetry breaking which de- 
cays exponentially versus the perpendicular distance z. 
The penetration depth ^ is set by the thickness of the 
layer of Mn dopants deposited at the surface; we take 
^ ~ Inm, assuming they penetrate several atomic lay- 
ers into the surface. In a real system the T-breaking 
term may vary depending on details of the system. Here 
we take a simple T-breaking term. The Hamiltonian 
H = Hq -\- Hi can still be written in the Dirac form 
H[0{z),k] = {e{k) - fi^ eV{z))I-hj:ada{k,0)r-. All 
the z dependence can be absorbed into the parame- 
ter by defining z{0) = -^ln|l - 0/7r\. da{KO) = 
d^a^ (k, 0) + ^^,3^0 sin[^] - Sa,5Mo{l + cos[0{z)]). FoT this 
model Hamiltonian, a^^'^^^ vanishes. / in Eq. (7) can 
be shown (see Supplementary Materials) to be given by 

/ = [_ d^kdOetahcddtdxdadydhdzdcdedd, (10) 

where the integral is carried over the outside of electron 
Fermi pocket, denoted by FFS and d = d/\d\. As dis- 
cussed earlier, can be viewed as the fourth momentum, 
playing the same role as the three-dimensional momenta 
kx,y,z, and the three-dimensional phase space Fermi sur- 
face can be defined in the fe, space. 

Numerical results of the surface Hall conductivity are 
shown in Fig. 2 as a function of the surface gate car- 
rier density, rig = CVg, for different strengths of the gap 
induced by the magnetization, mo. Note our theory is 
only valid when the surface is gapped. The Hall con- 
ductivity saturates to the quantized value /2h in the 
limit that Zd ~ Tig/niy ^ ^. Physically this is because 
the T-breaking, and thus the nonzero Hall conductivity, 
only occurs in the range of z < ^. For ^ <C the system 
is insulating in the range of 2; < ^, so that the surface 
Hall conductivity remains at the quantized value of TI 
surface. As expected, in the limit of low carrier den- 
sity the Hall conductivity also approaches e^/2/i, which 
is dramatically different from an ordinary semiconductor 
in which the intrinsic Hall conductivity always vanishes 
in the limit of zero carrier density. 

If observed, such a gate tunability of Hall conductivity 
in a hulk crystal provides a unique signature of doped TI 
compared to ordinary bulk semiconductors. We note that 
while disorder may add extrinsic contributions to the sur- 
face Hall conductance that we calculate here, these con- 
tributions vanish as the gate voltage is increased and the 
quantized e^/2/i value is recovered. Our work provides 
quantitative estimates for achieving the quantized value, 
and by comparing with measurements, can be used as a 
way to estimate the strength of extrinsic contributions 
to the surface Hall conductance. Following the approach 
of Ref. [1] for topological insulators, our results can also 
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be extended to general dimensions. Such generalizations 
lead to a framework to study topological phenomena in 
Fermi liquids in generic dimensions, which we investigate 
in future work. 
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I 

SUPPLEMENTARY ONLINE MATERIAL 

The purpose of these Supplementary Materials is to prove the following results (eqns. (7) and (8) of the main text): 

f = f d^kd^Trir^^'' r^^~' r^^~' r^^~' r^^~' 1 

•' ~ 15 J (27r)5 dq'' dqP dq" dq- ^ 

= ^ J^/^^ '^^'^ (J'abTcd) e''""', (11) 

and 

/niodl = ^^^^ CcsiA'''}. (12) 

The quantities A and are defined in the text and will also be defined below, q = (co',k) is a five- vector including 
the four-component spatial momentum k and frequency cj, and G = +ikgn(^ ) non-interacting single- 

particle Green's function; Pa is a projector onto band a and = — where fi is the chemical potential. For a 
metal, / is not quantized; despite appearances, it is not a second Chern number because the definition of the matrix 
T is such that its dimension changes at diff'erent points in momentum space, depending on the location of the Fermi 
surface (s). 
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The sections below are organized as follows. We begin by first manipulating the expression in terms of Green's 
functions into a more useful form. Next, we point out that even though / is not invariant under arbitrary deformations 
of the band structure, it is invariant under a certain class of deformations that we discuss. Following this, we prove 
some mathematical identities that the Berry connections of the Bloch bands must generally satisfy. Finally, after 
developing these ideas, we show how to use them to prove eqns. (11) and (12). We end with a discussion of the 
calculation of / in the simpler case of the Dirac models (Eq. (9) of the main text), which have additional symmetries. 



PRELIMINARY MANIPULATIONS OF / 

We have 

G"^(k,cj) =u;-iJk + ^^sgn(iJk), 
d^G-^ = 1 

d,G-' = d,H^ ^ {a\diH^\p) = diU^)So^P + (^^3 - ^a){a\di\p) 

G(k,u;) = y -—, (13) 

^ ^ ^a;-e, + z^sgn(e,)' ^ ^ 

where di = ^ is a derivative with respect to a momentum coordinate. Inserting a complete set of states for each k, 
we have: 

^ = (2^^'''' / ^ AfA^-'Al'At'^dldpd.ds, (14) 

where we have set A^^ = {a\diG~^\^) and = ^ — + iSsgn{^a)' 

Now observe that the above sum vanishes whenever at least any three of a, /3, 7, S are equal. For example if 
/3 = 7 = ^, then the summand is dP^(fi^A^^ A^-^ A^^ A^^ ^ which vanishes upon contraction with the epsilon tensor. 
Similarly, if a = /3 = 7, the summand is d'^dsAf^A^^A^^Af^, which also vanishes upon contraction with the epsilon 
tensor. 

Thus we are left with 3 types of terms: 

l^e'^"^ J d'kd^ ( E + E + E ) AfA^-'Al'At'^dld.d.ds (15) 



3 (27r) 



In X]^^^^, the sum is over states (a, /3, 7, (5) that are all different. In X]a^7(55 ^^^^ states where exactly two 

of the indices are equal and the others are all different. Finally, Y^'^^^^ is over states where two pairs of indices take 
the same value (e.g. a = /3 and 7 = ^ but a 7^ 7). 

Let us first show that the last sum gives zero contribution. Its contribution is proportional to: 



/ 



It is easy to see that the second term vanishes because of the epsilon tensor, and that the first and third terms 
cancel each other after indices are relabelled (to see this, note that e^*-^^ = — e*-^^^, since cyclic permutations have odd 
signature in even dimensions). 

Therefore, we see that the third sum in (15) vanishes. Now consider the second sum in (15), which we will call B\ 



n 

^ e'^^^ J d^kdujAfA^^Al^Af^dld^d^ds (17) 



^ ^ a/3'y5 

Recall that is a sum over states where exactly two of the indices are equal and the rest are all different. There 

are exactly 6 ways this can happen; let us write each one out: 

B = Ba = p + Ba = ry + Ba = S + B^=ry ^ B^^S ^ ^7 = 5, (iS) 
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where the subscripts indicate which terms in the sum are to be set equal to each other. We find that Bo,=^-\-Bo,=s = 0- 
We also find 



B^=j + Bj=s 



iiidicates a sum over states where a, /3, 7 are all different, c.c. denotes complex conjugation. Note that 
{Af y = -Af ^ if a ^(3 while Af^ = di^a is real. 

Therefore, we see that / is given by the following sum of three terms: 

cx.(3'y5 

y{k) = e'^>^^J2^l^^'sAfA^"AfA'r f dujdldpds{d^ - d^). 

t 

^(k) = -ie'^'^' dmUifiio^Mf + AfAfAt,'') J dLodldjd-y, (20) 

where we have set iA'^^ = {a;'k\di\P;'k) and £^a/3 = £,a ~ C^- The important point to note here is that x(k) and y(k) 
do not depend on the slope of the the dispersion, di^a, while z{h) does. 

Note that in the sum in z{k.), all three bands a, P and 7 must have different energies. Consequently, for a special 
class of models such as the Dirac-type models, which have only two bands with distinct energies (each of which 
generically has a degeneracy), z(k) must be zero. In these cases, it is easy to show that the integrand in the remaining 
terms are independent of the energy spectrum and are proportional to TrJ^ A J-". Thus a direct and explicit evaluation 
of x(k) + ^(k) reveals eq. (Fl). 



Topological invariance of / upon smooth deformations of G 

For general band structures with three or more distinct energy levels, the proof of eq. (11) is more nontrivial and it 
is essential to exploit the topological invariance of / to deformations that preserve the band structure near the Fermi 
surface. Under a smooth deformation of G to G -\- SG, the variation of / is given by 



^ ^ijkl 



3 (27r)5 



a/3S 

2 



J2 J d^kdioAfA^"AfAf"dldpds{d^ - d0). 



i_gy/c; ^ J d'^kdujiAfAfA^-'A]'' - c.c.)dldld^ (19) 



3 (27r) 



Sf = —e^-P-- j -^d^Tr[{G-^5G){d,G-^G){d,G-^G){d,G-^G){d.G-^G)\. (21) 

Note that a deformation of the band structure that changes the location of the Fermi surface is not a smooth 
deformation, because it can move the poles of G from the upper half plane to the lower half plane (or vice versa), and 
these are not small changes. Therefore, a small deformation of the band structure corresponds to a small deformation 
of G provided that it is away from the Fermi surface. Thus we can only deform the band structure if we keep a small 
window near the Fermi surface unchanged. Observe that since the integrand in (21) is a total derivative, it will vanish 
if 5G vanishes at the boundaries of the integral. Therefore, we can deform an arbitrary band structure (Fig. 3a) into 
a fiat band dispersion (Fig. 3b), where the band structure is unchanged in a small window near the Fermi surface. 
We may take the limit that the size of the window goes to zero at the end of the calculation. 



8 



(a) 



(b) 





FIG. 3: (a) Cross-section of a generic band structure, with the Fermi level lying within a band, (b) Deformation of band 
structure in such a way that / is left invariant. The filled and empty bands are deformed to a flat band, while the partially 
filled bands are deformed to flat bands everywhere except a small window near the Fermi surface, where they are left unchanged. 



FS stands for an arbitrary 3D Fermi surface in the 4D Brillouin zone. J^(^) is the curvature associated with M 
bands, including m bands below the Fermi surface labeled by 7 and M — m bands crossing the Fermi surface labeled 
by P. Our proof will require us to be able to set ^p(k) > ^^(k) for all k, so (22) is true only for those bands for which 
this is possible. 

First we observe that for an insulator, 



where is the field strength of the Berry connection associated with the filled bands. This formula is easy to prove 
in the case that the bands are all flat, with the empty/filled bands having energy or ^i?, respectively [1]. For a 
general band structure, it is argued that / is invariant under deformations of the band structure, and so its value is 
given by deforming the energy dispersions to the flat band model. This means that the terms in (20) that depend 
on energy must vanish, imposing various constraints on the properties of the Berry connections. (22) is one of these 
constraints; here we will derive it, along with some others as well. 

First we consider a band structure with flat energy dispersions ^p, and ^p, each of which may be multiply 
degenerate. When and are empty, it is easy to show that 



Identities 



For the proof of eqn. (11) and (12), here we will prove the following identity: 




(22) 




(23) 



a;(k) + y(k) = - 6 * 2iTi{AEFEF + ApppF + Apfef + Aefpf) 



oo 




(24) 



We have set 






(25) 



(q!,/3,7,5)G(A,B,C,P>) 



where the indices A, 5, C, D can refer to P, or F. When the bands are all flat, z(k) = 0, so / = — 
y). Since the second terms (whose coefficients depend on the energies) must vanish, it follows that 





(26) 



Now consider a general band structure for an arbitrary insulator that has M filled bands. We know that we 
can deform the bands however we like without changing /, as long as no band crosses the Fermi level during the 
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FIG. 4: (a) Cross-section of a generic band structure with M filled bands, (b) Deformed band structure: the empty bands are 
deformed to a flat band with energy A, m filled bands are deformed to flat bands with energy —A, and M — m filled bands 
are deformed so that the energy is flat everywhere except in a small window of width 2A where it varies linearly, as discussed 
in the main text. 



deformation. Let us deform them in the following way: consider a band structure where there are bands with energy 
(^E that are empty, m bands with energy that are filled, and M — m bands with energy that are filled. 
and will be flat bands, with energy ±A, respectively, will be flat everywhere except a small region in the BZ, 
where it varies linearly (See Fig. 4). More specifically, pick a boundaryless 3-manifold, called FS^ and suppose that 
varies linearly with distance from FS in some shell of thickness 2A around FS. In this case, Z = J^^ d^kz{]<.) is 
given by 

Z = - ^e^Jk^ E E / d'm^p^Pa^^p^aMf^k^^r + A';^AfAfn [ dujdldld,, 

p JfS±K J 

= -ie^jkiY^ I d^kduj^FP^EF^PEd^^pdldEdpidE-dp) E {AfA^^Ar^A';^AfA['') 
p Jfs±a i3eE,-feF 

=e'^^' E / d^kdujiFPiEFiPEd^ipdidEdF{dE - dp)]: E (^f^-^fcT + A^ :f^C) 

o JFS±A ^ 

d'^kdujV'di^P^pp^EF^PEdpdEdp{dE - dp) 



FS±A 

=27ri / d^kV%^p^PF^FE^EP[jT^ + - tAt- - TT^]. (27) 

JfS±A SFPSF£; SFPSFF '^PE'^PF <ipE<iPF 

where we have set F*(k) = ^^p^^^piA^^ J-'^f' + X-^J^^j^). For the sake of being explicit, we set the energy 
dispersions to be: 

= -^F = A, 

r A/2 if kG {FS±k) 

^p(k) = <^-A/4 if k e (BZ\(FFS' + A)) (28) 

[ -3A/4 if k e \fFS - A) 

Here we have set k to be the distance in the direction n normal to FS. 

^PF = A(l + Cp/A) = A(l/2 + K/4A) 

^EP = A(l - Cp/A) = A(3/2 - /t/4A) (29) 
Z = 2« A f d'kV'^^pE^EE^EPljT^ + tA^ - tA^ - jA-]- (30) 

JfS±A ^FP^FE ^EP^EF SPFSPF <;pE^PF 

After some simplification, this becomes 



'^dK'|^^d3fc||y"(AK',k||)(i + ^)(^-^)x 

1 1 1 



^4(1/2 + k'/4)2 4(3/2 -/t74)2 (1/2 + /tV4)2(3/2 - k'/4)2 (3/2 - k'/4)3(1/2 + /t'/4)^ 



(31) 
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where we have made a change of variables = k/A. Thus, it is clear that 



lim Z (X V • n 

A,A^O Jpc^ 
A/ A fixed 



/ V-h (32) 
Jfs 



Now observe that 



f=^J d^kTrF AF + ,5/°"* + 6r, (33) 

where Sf^^^ arises from the contribution of integrals outside of FS'zb A, while Sp^ arises from contributions of integrals 
over FS ± A: 



^j^out . '^^^ 1 



Q fo f]{-^Ep/^FET if + / 1 d^k[{ApEPF - ApEFE - Apfef)] 

^ \^^) ^2 \Jbz\{ffs+k) Jffs-aJ 

• 2 1 ^ r 

Sr = ^j^^STTi V(-1/^F^)" / d^k^^p(k)[ApEPF - ApEFE " Apeef] + Z (34) 

^ \^^) ^2 JfS±A 

For any choice of A, it must be the case that 

because / is actually given by the first term of (33). First observe that 

lim SP'' = lim Z. (36) 

A,A^O A,A^O 
A/ A fixed A/ A fixed 

This is because in that limit, it is clear that X]^2 /f5±a ^^^[(^^^(^)/^^^)^(^^^^^ ~ Apefe — Apfef)] 0, 
because the integrand is finite while the integration region becomes vanishingly small. Furthermore, 

lim Sr^' OC Ti-^Ep/M'' / d^k[{ApEPF - Apefe - Apfef)] = 0. (37) 

A,A^O J 
A/ A fixed ^^=2 

The last equality above follows from eqn. (26). From the above equations, it follows that 

lim Z (X [ V'h = (38) 

A,A^O Jpg ^ ^ 

A/ A fixed 

Now observe that 



/ y.noc6^^'^V/ 4''y^''^T^.^y^^'' = 0^ (39) 
Jfs ^ Jfs 

which proves the identity (22). 



Proof of Eqn. (11) 

Armed with the identity (22), here we will prove eqn. (11) by deforming an arbitrary band structure - in the case 
where the chemical potential lies inside a band - in such a way that keeps / unchanged while making it amenable to 
calculation. 

Consider the following deformation of the band structure, which keeps / unchanged because it does not deform 
the partially filled bands in a window near the Fermi surface (see Fig. 5). As in the previous section, we will use the 
variable k to denote the distance from the Fermi surface, and to be the directions parallel to the Fermi surface. 
Choose some small energy A near zero so that there exists Ai and A2 where ^p(a^ = Ai) = A and ^p(a^ = A2) = — A. 
Deform the empty and filled bands to have a fiat dispersion, with energy zbA. Then, deform the partially filled band 
for > Ai to have a flat dispersion with energy A. Finally deform the partially filled band for < A2 to have a flat 
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|Ai -A2| 



FIG. 5: (a) Cross section of generic band structure with partially filled bands, (b) Deformed band structure. The bands 
are deformed to flat bands with energy ±A relative to the Fermi level everywhere outside of a small window near the Fermi 
surface. Inside the small window near the Fermi surface, the bands are unchanged. The limit A ^ is taken at the end of the 
calculation. 



dispersion with energy —A. For A2 < < Ai, the dispersion of the partially filled band is unchanged. The region 
where A2 < < Ai will be denoted as FS ± A. 

Now we can split up the integrals over momentum space into a sum of integrals: 

ZTT^ 1 

/ = 3~ (27r)^ [{^out + yout + Zout) + {^in + + ^in)], (40) 

where Xout = Ibz\{fs±a) Xin = Jjrs±A and analogously for Y and Z. Zout = because it depends 

on 9^Q,, and 



7r2 1 



3 (27r)5 

(41) 

where J^^^) is the Berry curvature for the M filled bands outside of FS^ while is the Berry curvature for the 

M filled bands and the P partially filled bands whose Fermi surface FS we are considering. 

Observe that in the limit A ^ 0, Xin + Yin because the size of the region of integration goes to zero while the 
integrand stays finite. In the following, we will show that 



lim Z (X e'^ 



^jky^f ^(M);P,^(M);,P^Q^ ^^^^ 
17^ JFS 



IFS 

which will show that in the case of a single Fermi surface, 

Jbz\ffs Jffs 

The generalization to multiple non-intersecting Fermi surfaces is straightforward. 
Now consider Z: 

P ^7 JfS±A J 

= _i,^jkiY^ I d''kduJ^FP^EF^PE^^^pdUEdF{dE - dF) ^ {Af A^^ A^ ^ A^^^ Af Af'^) 

=e^^^' V / d^kduj£^Fp£.EF(,PEd,(,pdldEdF{dE - ^^f) J Y^^A^^ + A^ :F^,^) 

P Jfs±a ^ ^ 

= / d^kdujV'di^p^Fp£,EF^PEdpdEdF{dE - dp) 
Jfs±a 

=2m [ d^'kV'd.^p^PF^FE^EPl-^. 

Jfs±a sfpv 



^FP^FE ^EP^EF ipE^PF (,PE^PF 

2np(-4^--^^)], (44) 

t,PEt,pE <;PE^PF 
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We have: 

^PP = A{1 + X) 
^EF = 2A 
^BP = A(l-X), 
X = Cp/A. (45) 

Thus 

Zoc^ [ d^kV-d^pI{k,A,), (46) 
^ Jfs±a 

where 

/(k, Ai) = (1 + X)(l - x)[— ^— + —r\^ - r. — ^Z. , + 7^ — , + 2np( ^ 



(1-X)(1+X)3 

Simphfying: 



4(1+X)2 4(1 -X)2 (1-X)2(1+X)2 (1-X)(1+X)3 "(1-X)3(1+X) 
^ )]• (47) 



ii-^ if X(k)>0 



/(k,A.)^|aggi^, if x(k)<o 

In the hmit A ^ 0, which by definition also takes ^ while keeping A/A^ fixed, we see that the size of the 
integration region is proprotional to Ai — A2, while /(k, Ai) is non-singular. It follows that 

lim Z (X [ V-d^poc [ V'h, (49) 
A^o J^s Jfs 

where the latter proportionality follows because has non-zero slope at the Fermi surface only in the direction 
normal to the Fermi surface. From the definition of V, (42) follows. 



PROOF OF EQN. (12) 

In the last section, we proved the highly non- trivial relation 

/ = 3^ y d^kTr (TabJ'cd) 6"''='. (50) 

The dependence on the chemical potential is implicit in the definition of the field strength ^(k). J-(k) is defined 
to be the field strength of the Berry connection for the bands that are filled at the point k in the four-dimensional 
Brillouin Zone. Depending on the location of the Fermi surface(s), the number of filled bands will be different at 
different k-points. This is not the typical formula for a second Chern invariant. 
Here we prove a more mathematical result: 



(51) 



where A^^'^ is the Berry connection involving only the partially filled bands that are part of the ith Fermi surface 
FSi. 

In order to prove this, let us first concentrate on the case of a single Fermi surface. We may rewrite / as 

^=3iL^„.^(^-°'-"-)'-'"+3iL^(^S'^iJ')'-". (52) 

where FFS indicates the filled Fermi sea, ie the region in momentum space inside the Fermi surface, while BZ\FFS 
indicates the region outside of the filled Fermi sea. J^^^^ involves the completely filled bands, while J^^^"^ involves the 
completely filled bands and the partially filled bands. 
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Now observe that FahFcd^^ is locally a total derivative: 

Fabi^cde"'"^' = ^dMkdcAa + |^,A,Ad)e'"'^^ (53) 

which implies: 

/ mod 1 = ^ Jj£cs{A^'^) - (54) 

where A^^"^ is the non-Abelian Berry connection for the completely filled bands, A^^"^ is the non-Abelian Berry connec- 
tion for the completely and partially filled bands, and C^^{A) = e^^^TT[A^duAx + ^A^A^^Aa] is the Chern-Simons 
form. Since the Chern-Simons form is non-linear, this subtraction has a cross-term: 



/ mod 1 = 




The second term vanishes by Eq. (22), so that we obtain 

/modl = -i^ / CcsiA^''^)- (56) 

^TT^ Jfs 

The generalization to multiple Fermi surfaces is straightforward and yields (51). 



/ FOR DIRAC -TYPE MODELS 

Consider the Hamiltonian 

H{k) =e{k)I ^daT% (57) 

where F^ for a = 1, 5 are the 4 x 4 F matrices; their relation to the physical properties of Bi2Se3 are described in 
[2]. The eigenvalues of this Hamiltonian are ^±(k) = e(k) zb |c?(k)|; each band is doubly degenerate. P± = ^{l±d-T) 
projects onto the bands with energies The Green's function is: 

P 



- ^+ + z^sgn(^+) - ^_ + i^sgn(^_) 
^ - (e+ + e-)/2 + iSjsgniM + sgn(e-)]/2 + d • F[(^+ - e_)/2 + zS{sgn{^_) - sgn(e+))/2] 
(^ - C+ + ^^sgn(C+))[u; - + iSsgn{^-)] 
uj-\-d-r 



(^ - ^+ + ^^sgn(^+))(a; - ^- + i(5sgn(^_)) ' 



(58) 



where a; = cj - (^+ + ^_)/2 + zJ[sgn(^+) + sgn(^_)]/2 and d = d[{^^ - ^_)/2 + i^(sgn(^_) - sgn(^+))/2]. Substituting 
into the formula gives 

_7r2 rd^kdw 1 

/ - - y^ii/c/ J (2,r)5 [{c, - ^+(k) + iSsgn{^^)){u; - e_(k) + iSsgn{^-)W 

X TrliCb + d-Tfdi{e{k)I + daT''){u; + d-T)dj{e{k)I + daT''){u; + d-T)dk{e{k)I + ^iar")(c:^ + d-T)di{e{k)I + 6/,F^)] 

(59) 

Observe that any term with two or more factors of de wih clearly vanish under the epsilon tensor. The term with one 
factor of de can be shown to vanish. The term with no factors of de gives: 

_ TT^ f d^kduj didadjdbdkdcdidd 
/ - y eiiw J -j^ 



)5 [(^ _ + iSsgn{C+)){uj - + ifcgn(e-))]5 

X Tr[{oj + J - r)2r«(w + d- r)r'')(w + J • r)aftr<=(w + J - nr"*] (60) 



14 



A (eV A) 


3.33 


Bo (eV A) 


2.26 


Co (eV) 


-0.0083 


Ci (eV A^) 


5.74 


C2 (eV A2) 


30.4 


Mo (eV) 


-0.28 


Ml (eV A^) 


6.86 


M2 (eV A^) 


44.5 



TABLE I: Parameters of the model Hamiltonian for Bi2Se3. 



The trace simphfies to —4:etabcdd^{(^^ — d?)^ . Thus: 

-2 d'^kduj .......... 1 



f =^-^^ijki I -j7z—rr^tahcd(fdidadjdhdkdcdiddY. — —tt. xw r — — r xxio (61) 



If the chemical potential lies in the conduction band, this is: 



_7A f d'^k d^djdgdjdbdkdcdidd 

J — A ^ijkl I \4 ^tabcd 



4 ''''JszXFFsi^^r'^"'^^ \d\' 
3 
8^ 



lBZ\FFS 

If the chemical potential lies in the valence band, we have: 



/ d'^ketabcdd^dxdadydbdzdcdedd (62) 

Jbz\ffs 



IFFS 

For Bi2Se3 with a surface T-breaking field, we have 



f = -^[ d^ketabcdd^da^dadydbd^dcdujdd (63) 
Stt^ Jffs 



e(k) 




di(k) 


= Aoky 


rf2(k) 


= -AA;^; 


c;3(k) 


= mo sin ^ 


6^4 (k) 


= Sofc^ 


c?5(k) 


= -Mo cos 61 + Ml A;^ + M2A;j| 



(64) 

where we have only kept terms up to quadratic order in k. The parameters are given in Ref. [2]; we list them in Table 
I for completeness. 



EFFECT OF TOP GATING AT THE SURFACE OF A DOPED TOPOLOGICAL INSULATOR 

Here we would like to develop a calculation of the effect of gating at the surface of a doped topological insulator. 

We assume that the chemical potential is slightly above the bottom of the conduction band, and that the conduction 
band bends upward so that at the surface, the chemical potential lies in the bulk band gap. We also assume that a 
gap m has opened in the Dirac surface states and that the chemical potential is always inside of this gap. 

In the simplest depletion length approximation, we assume that within a depletion length Zd from the surface, there 
are no carriers; if the bulk doping charge density is — en5, where rib is the bulk carrier density, this means that there 
is a constant positive charge density erib in the depletion region. In order to maintain overall charge neutrality, there 
must be a compensating negative charge near the surface. Let us assume that all of this compensating negative charge 
comes from the application of a top gate. In practice this is not entirely true, but it may be a good approximation. 
With such an approximation, we have: 



VgC = eribZd, 



(65) 
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where C is the capacitance per unit area of the region between the metahic top gate and the surface of the material. 

This negative charge is accumulated in a metal gate a distance d from the surface of the material. Thus, the profile 
of the charge density is: 

p{z) = —CVg8{z ^ d) ^ en})Z^ z < Zd (66) 

Note that 2: = is the surface of the material, and z = —d is the position of the top gate. The potential V{z) is 
therefore given by 

cfv/dz^ = -p(z)/ieiz)eo), (67) 

where e(2:) = for — < ^ < is the dielectric constant of the dielectric that is in between the metal gate and the 
surface of the material, and e{z) = e for z > is the dielectric constant of the topological insulator. 

We are interested in V{z) for z > 0. For simplicity, let us assume that d ^ 0. The electric field in the region 
< z < Zd is: 

The potential is then 

V{z) = - r E{z) = -^zd{l - -f = - "-^r (69) 

The potential energy associated with this is eV{z). Note that Zd = VgC/enij = Ug/ni), where Ug is the number density 
of carriers on the top gate. 

Thus, we would like to fix the penetration depth of the magnetization, ^, the gap of the Dirac surface states m 
(which is induced by the magnetization), and the bulk density 715, and then study the surface transverse conductance 
as a function of the gate voltage Vg. For Bi2Se3, we take e ^ 100, so eeo ~ 8.8510~^^F/cm. 
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